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THE  CALCULATION  OF  eAt  WITH  SOME  APPLICATIONS 

1 . Introduction 

Throughout  this  paper  A will  be  an  n x n matrix  with  real  or  complex 

elements,  having  f(x)  as  its  characteristic  function  and  eigenvalues 

<*1 , cij,  ...  a (not  necessarily  distinct).  In  this  paper  we  calculate 

At 

the  exponential  matrix  e , specify  some  properties  of  certain  matrices 

At 

that  must  be  determined  in  order  to  describe  e and  finally  indicate 

some  applications  of  this  calculation. 

At 

2.  An  Initial  Value  Problem  and  e 

In  [3]  eAt  is  obtained  as  the  unique  solution  to  the  following 
initial  value  problem.  With  f(x)  as  specified  above  and  D e we  wish 
to  obtain  the  solution  to: 

f(D)G(t)  = 0,  G(t)  an  n x n matrix  with  elements  functions  of  t and 

(n-1)  n i 

such  that  G(o)  = I,  G' (o)  = A,  ...  G(o)  = An. 

At 

If  e is  defined  by  the  equation 
At  r AV 

e ' L ~ • 

n=0 

then 

eAt|  • I,  0(eAt) | = A,  ...  on-V*)!  ■ A"-1  , 

t=0  t=0  t=0 


and  f(D)eAt  = f(A)eAt  = 0 by  the  Cayley-Hami 1 ton  Theorem.  Therefore, 
At 

e is  the  unique  solution  to  this  initial  value  problem.  Suppose 
a-j , <*2 , ...  are  the  distinct  eigenvalues  of  A with  multiplicities 
u-| , ^2*  •••  us»  we  may  write 


(1) 


At 


■J, 


(C 


kl 


Ck2t 


C.  t 
kuk 


uk 


-1 


)e 


“k1 


3 


where  the  are  n x n matrices.  From  the  initial  conditions  we  have: 


1 = I ckl 

k=1  kl 


A ‘ (otkCkl  + Ck2J 


= kJ1  ^ak  1(:kl  + (n'1)ak  Ck2 


(n-l)!akn’yk 
(n-wk) ! (uk*l } ! Ckuk^  ’ 


from  which  the  n Ckj  can  be  determined.  We  note  from  the  system  (2)  that 
a solution  for  a Ckj.  will  be  a linear  combination  of  the  left  side  of  (2) 

and  therefore,  any  Ckj  will  be  a polynomial  in  A of  degree  at  most  (n-1). 

Being  polynomials  in  A,  the  commute. 

3.  Properties  of  the  when  all  roots  of  f(x)  are  distinct. 

We  will  not  use  the  system  (2)  to  completely  solve  for  the  n Ck^ , 

At 

rather  will  we  obtain  another  representation  for  e satisfying  the  initial 

£ akt 

value  problem  and  then  equate  coefficients  of  like  terms  "t  e ".  Before 
proceeding  to  this  representation  let  us  consider  the  simple  case  in  which 
all  roots  of  f(x)  are  distinct.  This  case  will  provide  insights  on  how 
to  handle  the  more  general  case  of  multiple  roots. 

If  all  roots  of  f(x)  are  distinct,  then  (1)  becomes: 

At  ^ t 

(l)1  eAt  = l C e k . 
k=l  kl 

The  inverse  of  (1 ) ' is 


-At  2 

e = l C.  ,e  k 
k-1  kl 


L 


Multiplying  (1)'  and  (1)"  and  using  the  conmutati ve  property  for  the  Ckl 
we  obtain 

n 9 n n (a.-a.)t 

(3)  1 - I Ckl  = l l C Ci]e  3 . 

k=1  Kl  k=l  j-1  kl  Jl 

k^j 

Equation  (3)  is  true  for  all  t,  however,  the  left  side  of  (3)  is  indepen- 
dent of  t which  suggests  that  Ck1  = 0 for  k ^ j,  i.e.,  the  Ckl  are 
pairwise  orthogonal.  (This  will  be  shown  below.) 

Applying  the  initial  conditions  to  (1)'  we  obtain: 

1 = cn  + c2i  + ...  + cnl 
A = + a2C21  + ...  + anCnl 


fln-l  _ n-lr  , n-1-  , . n-1- 
A ^i  Cn  + a2  C21  + ...  an  Cp1 . 

Since  the  c*k  are  distinct,  the  coefficient  matrix  for  the  system  (2)'  is 
a Vandermonde  matrix  and  is  non-singular.  Therefore,  the  system  (2)’  has 
unique  solutions  for  the  Ckl . If  we  solve  the  system  (2)'  for  , then 
the  coefficient  of  An"^  in  this  solution  is: 


(4)  (-1)' 


...  1 1 


Vl  Vl 


n-2  n-2 

ak-l  Vl 


/ al  a2 


n-1  n-1 


Both  numerator  and  denominator  of  this  coefficient  are  non-zero  Vander- 
monde determinants  and,  therefore,  C.,  (k  = 1,  2,  ...  n)  is  a polynomial 


. 


1 


in  A of  degree  precisely  n-1 . Moreover,  the  coefficient  of  A0"1  in  this 
solution  for  Ckl  is 

n 

1/  n (a,  - a.)  , 

j=l  * J 
j^k 

which  agrees  with  the  expansion  of  (4). 

Next,  suppose  we  wish  to  eliminate  Ckl  from  the  second  through  the 
n^  equation  in  (2)',  which  yields  n-1  equations  in  the  n-1  unknown 

matrices  , C21 » ...  Ck  -j  Ck+1  ...  C^.  We  do  this  in  order: 

’ 2 

subtract  ak  times  the  first  equation  from  the  second,  ak  times  the  first 

1 A.  L_ 

from  the  third,  to  finally  ak”  times  the  first  from  the  n . Deleting 
the  first  of  this  new  system  of  equations  we  obtain  n-1  equations  with 
the  unknown  missing.  The  left  members  of  this  new  system  will  be 
A - akI,  A2  - a2 1 , ....  An-1  - ak""*I,  each  of  which  has  a factor  A - a.1, 
as  does  any  linear  combination  of  these  left  members.  Therefore,  each 
of  the  solutions  for  C^,  C,^  , ...  Ck_-|  1 Ck+1  ^ ...  will  have  a 
factor  A - akI.  From  this  we  conclude  further  that: 
n 

(5)  Cji  = a^  n (A  - c^I)  , j s 1 , 2,  ...  n 

for  some  scalar  a,.  We  note  from  (5)  that  C.,  is  precisely  a polynomial 

J J 1 

of  degree  (n-1)  in  A with  leading  coefficient  a^  which  must  be 

(6)  a.  = 1/  n (a.  - a ) (from  (4))  . 

J £=1  J 1 

As  an  example  using  (5)  and  (6),  let  A be  any  3x3  matrix  with 
distinct  eigenvalues  , a^,  then; 


6 


At  (A  - ot2I)(A  - a3I)eait  (A  - aiI)(A  - ct3I)ea2t  ^ (A  - aiI)(A  - a2I)el3t 

(ct-|  - a2)(a-|  - a^)  (a2  - <*■] ) (<*2  - *3)  («3  - a-| ) (»3  - ci 2) 

3.1  Orthoqonal itv  and  Idempotencv  of  the  -CkT 

A further  conclusion  to  be  made  at  this  point  occurs  when,  using  (5), 
we  multiply  and  C.^  (i  f j);  this  yields 

cjicii  ■ Vt  j,  <A  - “<■>  j, (fl  - V> 

' ifj  ifi 

= a. a.  f ( A)  n (A  - a I)  = 0 
J 1 *=1  1 
£?i , if  j 

by  the  Cayley-Hami 1 ton  Theorem.  Therefore,  as  conjectured  earlier,  the 

C,,  1 are  pairwise  orthogonal  when  the  eigenvalues  of  A are  distinct. 

Using  the  first  of  the  initial  conditions  in  (2)'  and  this  orthogon- 

n 2 

ality  we  have  C .,  = C.,-,  ) C, , = C , j = 1 , 2,  ...  n,  i .e. , each  C ., 

J • J 1 ^1  ■ 1 J 1 J1 

is  idempotent. 

We  summarize  section  (3)  in  the  Theorem  I.  Given  A and  f ( X ) with 
distinct  eigenvalues  a^,  a2>  ...  ap  we  have 


At  - 

a)  eAt  = l Ck1e  k 
k=l  Kl 


b>  Ckl 


n (A  - a I) 

n 7 — r 

4*1  (c*k  ■ V 
ifk 


k = 1 , 2,  . . . n 


4. 


[ C.,  if  i = j (Idempotent) 

= » 11 

0 if  i f j (Orthogonal ) 
d)  All  Ck-|  are  polynomials  of  degree  n - 1 in  A and  they  commute. 
Properties  of  the  C.  . when  f(x)  has  multiple  roots. 
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•j 


f 

a 


4.1  Minimal  Polynomial  for  a given  matrix. 
First  consider  the  example: 


10  11 
0 10  0 
0 0 10 
0 0 0 1 


for  which  f (x)  = (x  - 1 ) . 


According  to  (1) 

eAt  = (C11  + C12t  + C13t2  + C14t3^et- 
Using  the  initial  conditions  we  obtain: 

cn  = Is  ci2  = (A  ' z)>  ci3  = IT  (A  * I)! 


C14  = 3T  (A  ‘ 1)3  ’ 


However: 


C13  = (A  * ^ = 


0 0 11 
0 0 0 0 
0 0 0 0 
0 0 0 0 


and,  therefore,  C-^  = = 0.  For  the  given  matrix  (x  - 1)  = 0 is  the 

characteristic  equation  and  (A  - I)A  = 0.  But  it  is  also  true  for  this 

matrix  that  (A  - I)J  = (A  - 1)^  = 0 and  A - I / 0.  In  this  case  A not 

only  satisfies  its  characteristic  equation,  it  also  satisfies  the  equations 

(x  - 1)  =0  and  (x  - 1)  =0.  For  a general  square  matrix  A,  the  lowest 

degree  monic  (leading  coefficient  equal  to  1)  polynomial  that  A satisfies 

is  called  the  minimal  polynomial  for  A.  In  the  example  above  (x  -1)  is 

the  minimal  polynomial  for  the  given  A and  for  this  matrix  (D  - l)VAt  = 0. 
At 

Our  solution  for  e should  then  have  been  written 

eAt  = (Cn  + Cut)el. 

In  general,  if  <]/(x)  is  the  minimal  polynomial  for  a given  matrix  A 
and  the  degree  of  ip( X ) is  m(<  n),  then  eAt  satisfies  the  initial  value 


problem : 


8 


<^(D)eAt  = 0 


and 


Dk(eAt)t=0  - Ak 


k = 0,  1 , 2,  . . . m - 1 . 


4.2  A Redefining  of  the  C.  . 

K J 


Let  the  distinct  eigenvalues  of  A be  a,  , a^, 

s Ui, 


a$  with 


ijj(x)  = n (x  - a.  ) K the  minimal  polynomial  for  A (written  in  factored 
k=l  k s 

form),  >_  1 , and  £ yk  = m(<_  n)  the  degree  of  x ) . Following  the 
k=l  y. 

theory  presented  in  [2]  we  define  X. (x)  = ^(x)/(x  - a.  ) ; then  X.  (x) 

uk 

and  (x  - a^)  are  relatively  prime.  Therefore,  there  exist  polynomials 
PkU).  q k ( X ) (degree  of  pk(x)  < uk)  such  that: 

(7)  Pk(x)Xk(x)  + qk(x)  (X  - a/k  h 1,  k - 1,  2,  ...  s. 

Then  define: 

E k ( X ) = Pk(x)Xk(x)  and 


(8) 


Ek(A)  = Ek  = Pk(A)Xk(A) 


k = 1 , 2,  ...  s. 


We  note  from  (8)  that  for  k f i Ek(x)E,(x)  is  a polynomial  multiple 
of  'i)( x ) and,  therefore,  E^  * E^  = 0 for  k f i (i.e.  the  E^  are  pairwise 
orthogonal ) . 

From  (7)  we  form  the  product: 

s uk  s s 

IT  qk(\)(X  - a.  ) = ifr(x)  H q,.(x)  - T (1 

k=l  K k k=l  k k=l 

Replacing  x by  A in  this  last  equation  we  obtain 


Ek(A)). 


s s s 

'{'(A)  n q.  (A)  = 0 = n (I  - E.  ) = I - l E.  , 
k-1  K k=l  K k=l  k 

which  follows  from  the  definition  of  Mx)  and  the  orthogonality  of  the 


Ek>  From  this  we  have 
(9) 


1 = l Ek- 

k=l  K 


9 


Multiplying  through  (9)  by  E , using  the  orthogonality,  we  have 


E = e2 

l i 


i = 1 , 2,  . . . S 


i.e.,  the  are  idempotent. 
Next,  define: 

V*)  = (x  - ctk)E(x) 


(10) 


and 


N,. (A)  = N.  = (A  - a.  I ) E.  (for 


1). 


If  uk  = 1 , then  is  a simple  root  of  ^(x)  and  for  such  roots 

N^(A)  = ( A)  = 0.  We  note  from  (10)  that  = (A  - a^I)  kEk  = xp( A) 

u,  -1 

and  f 0.  The  N^  are  said  to  be  nilpotent  of  index 

Additional  conclusions  from  definitions  (8)  and  (10)  are: 

(a)  All  E^,  are  polynomials  in  A and,  therefore,  commute. 

(b>  V'k ' V Vj  ■ Vj  ■ 0 <k  * 

(c)  We  have  the  identity: 

Ol)  A = f Ek(okI+Nk) 

k-1 

which  can  be  seen  as  follows: 

^ Ek(ak!  + V = J1  (\Ek + AEk  ■ \Ek) 


= y AE.  » A l E.  = AI  = A. 
k-1  K k=l  k 

At 

If  we  replace  A in  e by  the  identity  (11),  we  have: 
5 S S 


At  - a k=l 
e = e 


(J,  Ek(akI  * Nk))t 


(12) 


= e 


s 

y 

k=l 


LkV 


N, 


k-1 


S akt 


= I e {E,  + N.t  + 


(N.t)‘ 


(N.t) 


V1 


k-1 


21 


* P 


“k 


and  this  must  be  identical  with  (1)  i.e. 


y e {Ck1  + Ck2t  + ...  + C^t 

10 


k-1 


In  (13)  we  observe  that  we  have  only  the  Ek  (k  = 1 , 2,  . ..  s)  to  determine. 
These  can  be  calculated  from  the  definition  (8)  or  from  (2)  calculating 
only  the  Ckl  = Ek(k  = 1 , 2,  . . . s) . 

4.3  Summary  of  section  4 and  examples. 

We  summarize  this  section  in: 

s uk 

Theorem  2.  Given  A with  minimal  polynomial  d( x)  = n (x  - a,  ) , m < n, 

k=l  K 

we  have: 


At  ^ 

e = V E.e  {I  + (A  - ct.  I ) t + 


uk_1  uy'1 
(A  - a.I)  t 

+ , , K A 

+ ~TTi 


in  which  the  E.  satisfy  the  following: 
s 


(a)  I = l E. 


(b)  E.E . = 


( 0 if  i t j 
^ Ek  if  k = j 


(c)  (A  - 3|J)Ek  = Nk  are  nilpotent  of  index  uk 

(d)  all  Ek  and  Nk  are  polynomials  in  A. 

Up  to  this  point  we  have  determined  for  any  A(3  x 3)  with  distinct 
eigenvalues.  Now  let  us  complete  this  calculation  for  any  3x3  matrix 
A.  To  this  end  we  have  the  following  case^  and  calculations: 

3 

(i)  All  eigenvalues  of  A are  equal  to  a,  and  p ( X ) = f(x)  = (x  - a)  . 

In  this  case:  ? 

At  t M 

e = e (E1  + N] t + -yj— ) in  which  E]  = I and  * (A  - al). 


2 

(ii)  Again  all  eigenvalues  equal  a,  but  i|>(x)  = (x  - a)  . In  this 
case  E-j  = I and  e ^ = eat(I  + (A  - al)t). 

(iii)  ^(x)  = (x  - a).  In  this  case  E-j  = I and  eAt  = eatI  a scalar 
matrix. 

(iv)  i(j(  X ) = f ( X ) = (x  - a-j ) 2 ( X - a2),  a-j  t a2> 

In  this  case: 


At  al^  > a2^ 

e = e (E-j  + N] t)  + e E2  . 

We  can  solve  for  E-j  and  E2  (N^  = (A  ~ a-jI)E-|)  using  the  initial 
conditions  (2)  or  definitions  (7)  and  (8).  In  view  of  (7)  we  have 
(repl acing  X by  A) : 

E]  + qn  (A)  (A  - c^I)2  = I. 


However,  we  know  that  E^  + E2  = I and,  therefore,  E2  = q^(A)(A  - a-jl)' 
which  means  we  obtain  both  E^  and  E2  simultaneously  by  using  (7)  and 
(8).  Accordingly  using  (7)  and  (8):  (ax  + b ) ( X - a2)  + c(x  - a-j)2  h 

which  must  hold  for  all  X and,  therefore,  we  have: 


, b = 


- (a2  - 2a-| ) 


(a2  - a-j  ) (a9  - a-, ) 


2~  » c = 


(a2  - a-j ) 


2 ■ 


Therefore : 

E-j  = j (A  - a2I)(A  + (a2  - 2ai)I) 

(a2  - a,) 


E2  = j (A  - ajl)2  and 

( a2  ” al ^ 

- (A  - a,  I ) 

N,  = j-  (A  - a2I)(A  - 2a-j  I + a2I)  . 

( a2  - a-j ) 


Accordingly: 


e 


At 


a-i  t apt 

e 1 (Ej  + N-j  t)  + e Ep. 
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(v)  The  last  case  is  that  in  which  f ( x ) = (x  - a^)  (x  - a^)  as  in  (iv) 
but  <|;  ( X ) = (X  - a-j)(X  - c^). 

In  this  case  and  a ^ are  simple  roots  of  <i(x)  and  therefore, 


At 


Cl^  t 


(*2 1 


(A  - a,  I)  + 

a,  - a0  2 a0  - a. 


(A  - a-j  I ) 


by  virtue  of  the  results  in  section  3. 

In  view  of  the  example  in  section  3 and  the  5 cases  above  we  have 
At 

obtained  e for  any  3x3  matrix. 

It  is  of  interest  to  note  that  for  any  A(n  x n ) in  which  either 

(a)  A has  n equal  eigenvalues  or  the  opposite  extreme  (b)  A has  n 

At 

distinct  eigenvalues  we  have  that  e can  be  written  immediately  as: 


(a)  e 


At 


.at  V (A  - gl)ktk 


n-1 

l 

k=0 


k! 


If  in  this  case  ^(x)  = (x  - a)m,  m < n,  then  the  summation  would 
extend  only  to  m - 1 since  (A  - al)m  = 0. 
a.  z n 


At  2 akl 
(b)  e = ) e 

k=l 


(A  - a I) 

n 7 A-y  . 

j=l  uk  * aj} 
j^k 


At 

5.  Other  representations  for  e . 

In  [4]  and  [6]  e^  is  obtained  by  use  of  the  Lagrange-Syl vester 

interpolation  polynomial.  By  using  the  eigenvalues  of  A as  the 

At 

interpolation  points  we  obtain  the  form  of  e in  equation  (12).  In 

[5]  and  [7]  representations  of  eAt  are  obtained  in  one  case  in  powers  of 

A and  in  others  in  powers  of  A - I (a.  - eigenvalues  of  A).  In  any 

case,  if  all  these  representations  were  given  in  powers  of  the  same 

(A  - 8^1)  then  they  would,  of  course,  all  be  the  same. 

At 

Another  representation  for  e , which  has  applications  to  solutions 
of  first  order  linear  systems  of  simultaneous  differential  equations 
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with  constant  coefficients,  is  obtained  as  follows.  Suppose  we  have 
given : 

(14)  x'(t)=Ax(t) 

A an  n x n matrix  with  constant  elements,  x(t)  an  n x 1 vector  function 
of  t.  Suppose  we  have  found  a fundamental  set  of  solutions  for  (14) 
namely  x-j(t),  x2(t),  •••  xp(t).  Then  define: 

(15)  X(t)  = (Xl(t),  x2(t).  ...  xn(t))  , 

which  is  an  n x n matrix  whose  columns  are  the  elements  of  the  fundamental 
set.  X(0)  is  nonsingular  and  we  define 

(16)  G(t)  = X(t)(x(0))_1; 
then 

G(t)  = eAt. 

This  equation  follows  from  differentiating  (16),  which  yields 
G'(t)  = (x-j(t),  x2(t),  ...  x;(t))(X(0))*1  . 

By  (14)  G'(t)  can  be  written: 

G'(t)  = ( Ax 1 ( t ) , Ax2(t),  ...  AXn(t))(X(0))_1 
= AX(t)(X(0))_1  = AG(t)  . 

From  this  last  equation  it  follows  that: 

G(k)(t)  = AkG(t)  k=0 , 1,  .... 

which  in  turn  yields: 

G(k^(0)  = Ak  k=0,  1,  ...  . 

Moreover, 

f (D)G(t)  = f(A)G(t)  = 0 
by  the  Cayl ey-Hami 1 ton  theorem. 

Therefore,  G(t)  satisfies  the  initial  value  problem  which  is  also  satisfied 
At 

by  e . By  uniqueness  of  such  solutions  we  conclude: 

G(t)  = X(t)(X(0))_1  = eAt. 
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6.  Some  applications  of  e . 

From  the  last  remarks  in  section  5,  and  since  (X(0))~^  is  nonsingular, 
X(t)(X(0))_1  = eAt  has  columns  that  are  linear  combinations  of  the 
columns  of  X(t)  and  form  another  fundamental  set  for  the  differential 
equation  (14): 

X'(t)  = AX(t), 

Let  us  use  this  fact  to  solve  the  system 


x'(t) 


4-5  3 

2-3  2 

-1  1 0 


x(t). 


The  given  matrix  has  characteristic  polynomial  f(x)  = if<(x)  = (x  - 1;  (x 


Using  case  (iv)  in  section  4 with  = 1 and  02  = -1  we  have: 


2 -2  l 


0 00/ 


l -i  1 


1 -1  1 + 0 0 0 t ; + e 


1 -1 


0 0 


et(2  + t)  - e*1  -et(2  + t)  + 2e_t  e^l  + t)  - e_t 


-e1  + 2e_t 


et(l  - t) 


The  columns  of  the  latter  matrix  constitute  a fundamental  set  for 
the  differential  equation  and  its  general  solution  is: 
x(t)  = eAtx(0)  , 

where  x(0)  forms  the  initial  conditions  for  x(t)  (given  at  t = 0). 

At 

As  a matter  of  fact,  having  obtained  e for  all  3x3  matrices  A we 
have  therefore  obtained  fundamental  sets  for  all  systems  of  differential 
equations  X'(t)  = Ax(t)  with  A a 3 x 3 matrix  of  constants  and  x(t)  a 


3 x 1 vector  function  of  t. 
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r 


At  -At 

Knowing  e and  e we  define 
cosh  At  = l/2(eAt  + e'At) 
sinh  At  = l/2(eAt  - e'At). 


and 


Equally  well  we  know  e^At  and  e~^At  (i  = v^T)  and  define: 


cos  At  . l/2(e1At  * e-,flt) 


s1„  At  ■ ^(e1At  - e""L). 

As  an  example,  we  indicate  the  expansion  of  cosh  At: 

2 2 

S IT.  t S 


-iAt, 


3 3 

r,r 

cosh  At  = l cosh  ct^t( E^  + — +...)+  I sinh  a^t(Nkt  + -jy — + 


with  each  of  these  two  sums  terminating  with  the  term  containing  either 
wk‘2  uk_1 

t or  t , depending  upon  whether  is  even  or  odd. 

In  [1]  T.  M.  Apostol  considers  the  system  of  differential  equations 
Y"(t)  = AY(t) 

and  writes  the  solution  in  terms  of  two  matrix  functions 


C(t) 


- t2kak  - t2k  + l.k 

In  T2kTT  ‘ S(t)  = (2k  + 1)!  • 


k=0 


C(t)  is  precisely  cosh  /K  t and  S(t)  is  (*^)_1sinh  /K  t provided  that  A 

is  defined  and  nonsingular.  Clearly  one  would  define  /A  to  be  that 
2 

matrix  3 such  that  B = A.  It  turns  out  that  B is  not  unique  (as  one 

would  suspect),  in  fact,  there  may  be  as  many  as  2n  matrices  B such  that 
2 

B = A.  However,  these  B-  matrices  may  be  calculated  as  follows.  If 
A is  similar  to  a diagonal  matrix  so  that 


then 


A - i Di ag  {&■.,  ot^,  ...  ct  } 1*  , 


/A  = T"1  Diag  {+a]/2,  + z\/2 , ...  + c^72}  T. 
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.) 


If  A is  not  similar  to  a diagonal  matrix  and  is  nonsingular,  then 
can  be  obtained  as  follows:  From  (11) 


then 


/K 


r . 1/2  F Nkvl/2 

i.  ± ak  Ek(I  + 


k=l 


r'k  1 /2 

If  (I  + — ) ' is  expanded  by  the  binomial  theorem  then  this  expansion 
ak 


uk~] 


would  terminate  with  the  term  since  is  nilpotent  of  index  ut<- 

Knowinq  how  to  compute  JK  in  some  cases  we  then  have  for  these 
cases  the  solutions  to 


(a)  Y " ( t ) + AY(t)  = 0 
given  by 

Y(t)  = (cos  Jh  t)Y-j  + (sin  vTT  t)Y„; 
or 


(b)  Y" ( t)  - AY(t)  = 0 
given  by 

Y(t)  = (cosh  /A  t)Y.j  + (sinh  /A  t)Y^. 

In  (a)  and  (b) 

Y ( 0 ) = Yj  and  Y * (0)  = /A  Y£. 

As  a last  application  we  will  consider:  let  A be  a 3 x 3 matrix  with 

2 

characteristic  function  f(x)  = ij>(x)  = (\  - a1 ) (x  - a^)  (case  (iv)  in 
section  4),  and  suppose  we  are  qiven  the  nonhomoqeneous  svstem 
07)  x’(t)  = Ax(t)  + a^1  , 

in  which  x(t)  is  a 3 x 1 vector  function  of  t to  be  determined  and  a^  is 

-At 

a 3 x 1 constant  vector.  Multiplying  through  (17)  by  e yields 
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e’Atx'(t)  - Ae'Atx(t)  = D(e"Atx(t) ) = e'^a/1*. 

Integrating  this  last  equation  we  obtain 

(18)  x(t)  = eAt/te'ATa1eai  dx  + eAta2  , 

in  which  a2  is  a 3 x 1 vector  whose  elements  are  arbitrary  constants. 
From  case  (iv) 

a t ai t a2t 

(19)  eAt  = e (E-j  + N^t)  + e c Er 

Substituting  (19)  into  (18),  using  the  orthogonality,  idempotentcy  and 
nilpotentcy  of  E1 , E^  and  Ni  and  integrating  we  obtain 

a,  t N,  t E~  a,  t a~t 

x(t)  ■ e 1 { (E-j  + 4-)t  + + te  (E1  + l^t)  + e E2)a2 

as  the  complete  solution  to  (17). 
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